Control of beam propagation in optically written waveguides beyond the paraxial 

approximation 
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Beam propagation beyond the paraxial approximation is studied in an optically written waveguide 
structure. The waveguide structure that leads to diffractionless light propagation, is imprinted on 
a medium consisting of a five-level atomic vapor driven by an incoherent pump and two coherent 
spatially dependent control and plane-wave fields. We first study propagation in a single optically 
written waveguide, and find that the paraxial approximation does not provide an accurate descrip- 
tion of the probe propagation. We then employ coherent control fields such that two parallel and one 
tilted Gaussian beams produce a branched waveguide structure. The tilted beam allows selective 
steering of the probe beam into different branches of the waveguide structure. The transmission 
of the probe beam for a particular branch can be improved by changing the width of the titled 
Gaussian control beam as well as the intensity of the spatially dependent incoherent pump field. 



I. INTRODUCTION 

Diffractionless optical beam steering is of great inter- 
est in optical imaging, laser machining and optical com- 
munication. A gradient of the refractive index provides 
the simplest way to steer the optical beam. Such a re- 
fractive index gradient can be achieved by various phys- 
ical mechanisms such as thermal gradients acousto- 
optic interactions Q and electro-optic effects Simi- 
larly, light-matter interactions can be used to control the 
optical beam propagation through the atomic medium. 
More precisely, electromagnetically induced transparency 
(EIT) [| can provide beam-steering since refractive index 
changes significantly near the center of the transparency 
window [BlUl- 

Crucial obstacles in the above schemes are absorption 
and diffraction. The effect of diffraction is inevitable as 
an optical beam propagates through a medium. An op- 
tical beam can be considered as a superposition of plane 
waves with different wave vectors. Diffraction arises, 
as each plane wave acquires a different phase during 
its propagation. In free space, a light beam will typi- 
cally be distorted severely already after propagation over 
few Rayleigh lengths Q. For all-optical processing, this 
remains as a major obstacle in practical applications. 
In order to suppress or even eliminate diffraction, re- 
searchers have developed many proposals based on dif- 
ferent physical mechanisms, such as EIT [8l-fl5j. coher- 
ent population trapping (CPT) [H,[i3, coherent Raman 
processes [H, [l9[ , or saturated absorption [2(| ■ Most of 
the schemes employ suitable spatially-dependent struc- 
tures of the control field to prevent the optical beam 
from diffracting. Recently, it was found that alternatively 
Dicke narrowing induced by atomic thermal motion and 
velocity-changing collisions can be useful to eliminate the 
diffraction of a probe beam that carries an arbitrary im- 
age mm. 

A particular class of beam steering devices are so- 
called Y-branch waveguides, in which a single waveguide 



splits into two output ports |23H26l | . Y-branch waveg- 
uides can be used to divide a single beam into two sepa- 
rate branches with a certain intensity ratio. A desirable 
feature of such devices is a dynamical control of the light 
beam intensity at the different branches. 

Motivated by this, here we propose a scheme for all- 
optical beam steering in optically written waveguides. 
We facilitate a five-level medium driven by spatially de- 
pendent control fields, which was recently proposed as 
a method to achieve large refractive index modulations 
with low absorption [27j ■ We start by analyzing the light 
propagation through a single optically written waveg- 
uide. Comparing results obtained using a split-operator 
method assuming paraxial approximation to those from 
a finite-difference-time domain approach beyond paraxial 
approximation, we find that already in this simple case 
the paraxial approximation usually assumed in related 
studies does not provide an accurate description of the 
beam dynamics. Next, we consider an optically written 
branched waveguide structure consisting of two parallel 
beams crossed by a third tilted beam. This structure es- 
sentially consists of two coupled Y-branched waveguides. 
We find that the tilted beam can be used to switch the 
light propagation between cither of the two output ports 
formed by the two parallel light beams. Our numeri- 
cal results show that the coupling efficiency of the probe 
beam into the different branches can be controlled by 
changing the width and the angle of the tilted Gaussian 
control field, and the transmission of the output probe 
beam can be improved by increasing the amplitude of an 
additionally applied incoherent pump field. 

The article is organized as follows. In Sec. [Til we in- 
troduce the theoretical model, analytically derive the lin- 
ear response of the medium to the probe field from the 
master equation, and discuss the propagation equations 
and numerical methods to calculate within and beyond 
paraxial approximation. Sec. IIIII describes the waveg- 
uide structure written by a Gaussian control field inside 
the atomic medium, and discusses the possibility of con- 
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FIG. 1. The five- level type scheme considered in the anal- 
ysis. The probe field couples to transitions |1) <B> |4) and 
4} -o 1 5) and The two coherent control laser fields Q, s , il c are 
far-detuned from respective transition frequencies in order to 
split states |4) and |5) into suitable dressed states. An inco- 
herent pump control field is applied to |1) o |4) to achieve 
a population inversion. The control fields modify the probe 
field coupling to achieve high refractive index contrast with 
low absorption. 



trolling the optically written waveguide structure by an 
incoherent pump field. In Sec. IIV Al we compare nu- 
merical results for the light propagation through a sin- 
gle waveguide within or beyond paraxial approximation. 
Sec. llVBl presents our main results on the controlled light 
propagation in a branched waveguide structure. Sec. IVl 
summarizes the results. 



II. THEORETICAL MODEL 
A. Equations of motion 

We consider the five-level atomic system shown in 
Fig.[TJ This model has been proposed previously in a dif- 
ferent context for high index of refraction contrast with 
low absorption [27}. The two transitions |2) -f-s- |4) and 
1 3) <H- 1 5) arc driven by two far-detuned coherent laser 
fields with Rabi frequencies Q s and fl c respectively, and 
a weak probe is coupled to both transitions |1) <R- |4) 
and 1 4) |5) with Rabi frequencies f2 p i and f2 p2 respec- 
tively. An incoherent pump field is applied to transition 
|1) — > |4) to provide a population inversion. 

The electric fields of the three beams are defined as 

(r, t) = \e 3 (r, t) e, + c .c. , (1) 

where Ej (r , t) are the slowly varying envelopes and e*, the 
unit polarization vectors of the electric fields. The index 
j £ {p, s, c} labels the three fields, respectively In dipolc 
and rotating wave approximation, the Hamiltonian of the 



system in the interaction picture can be written as 

Hjjh = -(A pl - A s )|2)(2| - (A pl + A p2 - A s )|3)<3| 
- A pl |4)(4| - (A pl + A p2 )|5)(5| - (0 pl |4)(l| 
+ft p2 |5)(4|+ft s |4)(2|+ft c |5)(3|+H.c.) , (2) 

where A p i = oj p — w 4i , A p2 = oj p — uj 54 , A s = oj s - w 42 , 
and A c = uj c — W53 are the dctunings of the laser fields, 
and the Rabi frequencies of the fields arc defined as f2 p4 = 
Epfl^i ■ e p /2h, f2 p2 = E p jl54 ■ e p /2H, f2 s = E s p42 ■ e s /2h, 
and fl c = E c p^-e c /2h. Here, are the dipole moments 
of the respective transitions, and we have simplified the 
notation Ej(r,t) to Ej (j <G {p, s, c}). The master equa- 
tion of motion follows as 

p= l -[H I ,p]-C P (3) 

and Cp represents the incoherent contributions given by 
Cp = Cj lP + Cj 2 p + Cj 3 p + CJ 4 p + C d p + £Pp , (4a) 

q k p = (ii)oi p + p \m - mu\ p \s){k\) , (4b) 

£ p p = |(|l>(l|p + p|l)(l|-2|2)(l|p|l)(2|), (4d) 

where C~j k p describes spontaneous emission from \j) to 
|fc) with rate Tjk- The second term C d p models addi- 
tional pure dephasing for pjk with rate 7^, such that the 
total damping rate of this coherence is jjk = 7^. + (Tj + 
Tfe)/2, with Tj = ^ fe Tjk the total decay rate out of state 
\j). The third contribution CP p describes the incoherent 
pumping from |1) to |4) with rate p. 

The equations of motion for the density matrix ele- 
ments can be easily be derived to give 

P\\ = -ppu - i (QpiPii - & v iPi\) + r 4 ip 44 , (5a) 
p 22 = -iCl s (p 2 4 - p 42 ) + T 42 p 44 , (5b) 

piA = PP\\ - (r 4 i + T 42 ) p 44 - 1 ( — Slpiy0i4 - Sl s /5 24 

+SlplP4l + f2 s /3 42 + fl p2 p 4 5 — O p2 p5 4 ) + T5 4 p55 , 

(5c) 

P55 = —i (— — ^p2/945 + ^eP53 + Qp2P5i) 

- (r 5 3 + r 54 ) p 55 (5d) 

Pi2 = «(A S + *7 42 )p 42 + Kl s {p22 - pa) + ift P iPi2 

+ ift P 2P52 , (5e) 

P54 — — 754P54 - i (-^cP34 - ^p2P44 + ^plP51 + ^sP52 
-A p2 p 54 + f2 p2 p 55 ) , (5f) 

Pa = - (p/2 + 741) P41 - i (-^plPll - Q S P21 - A p ip 4i 
+ il p ip 44 - fi p2 /05l) , (5g) 

P34 = -743P34 - i (^plP31 + ^sP32 + (A c - A p2 )p 34 

+^p2P35 - ^cP54) , (5h) 
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P52 = -752P52 - i{-QcP32 - ^p2P42 - (A s + A p2 )p 5 2 

+f2 s/ 54 } , (5i) 

f>21 = - (p/2 + 721) P21 - i + A s )p 2 l + ^plP24 

-^sP4l} , (5j) 

P32 = -732P32 - i {(A c - A s - A p2 )p32 + ^sP34 

-O c p 52 } . (5k) 
The remaining equations follow from the constraints 

YsiPii = 1 and Pij = Pji- 

B. Steady state solution 

We assume the probe field to be weak enough to be 
treated as a perturbation to the system in linear order. 
The related zeroth [superscript (0)] and first-order [su- 
perscript (1)] contributions for are obtained as 
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(6a) 
(6b) 

, (6c) 
(6d) 
(6c) 

(6f) 



The detailed expressions for A, B are provided in Ap- 
pendix |XJ 



should be satisfied in Eq. ([7]) [23]. First, the popula- 
tion inversions on the two transitions should be matched. 
Second, the decoherence rates should be comparable. 
Third, the two-photon resonance should be fulfilled for 
the two transitions. These conditions can be expressed 
as r 44 (p 44 - pn) = r 54 (p 55 - p 44 ), 7 54 = 7 41 + p/2 and 
A p i + A p 2 = 0. A maximum refractive index modulation 
can then be realized if 741 = ±A p i — p/2. In this ideal 
case, the refractive index is increased by a factor of 2, 
with absorption canceled at the same time. However, it 
is hard to find a real atomic system in which these condi- 
tions are satisfied simultaneously. To relax the stringent 
conditions, it has been suggested in [23] to introduce the 
two far-detuned laser fields coupling |4) and |5) to auxil- 
iary states to induce Stark shifts. The Stark shifts of |4) 
and 1 5) together with modifications to the decoherence 
rates of the Stark sublevels can be controlled by tuning 
the intensities and detunings of the two laser fields. This 
way, the two transitions for the probe can be modified 
to one's advantage. This leads to the five-level scheme 
shown in Fig.[TJ in which the two upper states are coupled 
by two external far-detuned laser fields f2 s ,fi c . By ap- 
propriately choosing the parameters of these laser fields, 
one can achieve strong refractive index modulation with 
minimized absorption or even gain even in realistic level 
schemes. 



D. Propagation dynamics with and without 
paraxial approximation 

The propagation dynamics of the probe field is gov- 
erned by Maxwell's equations. Without applying the 
paraxial approximation, the wave equation for the probe 
field propagating along z direction can be derived as (28j 



dEp 
dz 



= (VI + ^)E P + -^x{u P )E p , (8) 



dz 



C. Interpretation of the level scheme 

To elucidate the role of the different couplings, we re- 
duce the level structure to a three-level ladder system, 
by switching off the two far-detuned laser fields. Then, 
the linear susceptibility for the probe simplifies to 



x(w P 



3N\l 
8tt 2 



r 4 i(p 44 — pa) r 54 (p 55 — p 44 ) 

) A p2 + «754 



i(7 



41 



(7) 



where N is the atomic density and X p is the wave length 
of the probe field. The terms in the right hand side of 
Eq. (UJ) describe the two individual contributions com- 
ing from the two two-level subsystems in the three- 
level ladder structure. To achieve the desired refractive 
index modulation without absorption, three conditions 



where k p = u> p / c with c the speed of light in vacuum and 
Vi = d 2 /dx 2 + d 2 /dy 2 . In Eq. ®, the diffraction of 
the probe, which is induced by the term icV 2 l E p /2u> p , 
may be reduced or even eliminated by the other term 
ik p x(uJ p )E p /2 if the real part of the susceptibility has a 
waveguide-like structure. 

To evaluate Eq. ([5]) without the paraxial approxima- 
tion, we make use of the finite-difference time domain 
(FDTD) technique [2^ provided by the software pack- 
age MEEP [3Cj . In order to compare the results to those 
obtained in paraxial approximation, we drop the term 
(ic/2ojp)d 2 E p /dz 2 which introduces effects beyond the 
paraxial approximation, and solve the resulting reduced 
equation using the Strang split operator method [311 ] . In 



both simulation techniques, all fields are chosen as con- 
tinuous wave fields. 
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III. REALIZING SPATIAL WAVEGUIDE-LIKE 
STRUCTURES 

A. Basic analysis 

From Eqs. (|6]), the linear susceptibility of the probe 
field can be written as 



3NXI 
8tt 2 



1 41P41 



1 54P 5 4 



„2 



(9) 
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We now assume a control field with spatial intensity de- 
pendence modeled as a Gaussian 



(10) 



where w s is the width of the intensity profile. This spa- 
tial dependence of the control field creates a modulation 
of the medium susceptibility in space, which can be used 
as a waveguide. In order to achieve maximum refrac- 
tive index modulation, the peak Rabi frequency of the 
laser field J7 S should meet the condition Sl s o = V2741 A s . 
Assuming these conditions, an example for the real and 
imaginary part of the susceptibility, which correspond to 
the medium dispersion and absorption, respectively, is 
shown in Fig. [2ja). Clearly, the real part of the suscepti- 
bility resembles a waveguide structure, whereas absorp- 
tion is low. 

The parameters in Fig. [2] are chosen to fulfill 
A p i, A s , 754, 741 > r 5 4,r 4 i,p. In this limiting case, 
Eq. ([6^) can be simplified to 



n 



pi 



-741 + «741 



Q'z 



(11) 



A s + 7 4 



In the central area of f2 s defined by Eq. (TU)) . the part 
f2g/(A s + 741) — 741 ~ 741 in the dominator of Eq. ([TT|) . 
and the refractive index is maximized as desired together 
with gain. At the same time, the transition |5) O |4) 
will give maximum refractive index together with absorp- 
tion. Those two transitions together result in a maximum 
refractive index for the probe with almost canceled ab- 
sorption in the central area of Q s as shown in Fig. [2ja). 
In contrast, in the two wings of Q s , the refractive in- 
dex is minimized for transition |4) -h- |1) with gain since 
f2g/(A s + 741) — 741 ~ —741. Simultaneously, a minimal 
refractive index with absorption is obtained for transition 
1 5) O |4). In total, a minimal refractive index with little 
absorption is obtained for the probe. 

In order to check the validity of the perturbation ap- 
proximation in the far-detuned case, we also obtained 
the probe field susceptibility to all orders by numerically 
solving Eq. ((5j) in steady state. For a weak probe field 
dpi ~ T41, it is virtually indistinguishable from the linear 
susceptibility shown in Fig. [5{a). 
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FIG. 2. (Color online) Real (blue solid line) and imaginary 
(red dashed line) part of the linear susceptibility obtained 
from Eq. <(9j as a function of the position in the presence 
of a control laser field f2 s with Gaussian intensity profile. 
Panel (a) shows results for a constant incoherent pump rate, 
whereas (b) shows results for spatially dependent pumping as 
explained in the main text. The inset in (b) depicts the de- 
tailed structure of the linear susceptibility in the central area 
which is similar to that in (a). Note that the profile of the real 
part in the central region forms a waveguide structure. Pa- 
rameters are chosen as: N = 1.4 x 10 27 m -3 , A p = 813.2 nm, 
d _ no r^rr„ ..d d = Q 3 GHZ] ^ = q 2 Q HZ] 

r 4 i = r 42 = 45 Hz, r 54 = 15 Hz, 

19.7 GHz, A s = 10.0 GHz, A c = 
f2 c = 5.2 GHz, w s = 1.0 fim, 



754 = 0.8 GHz 

„,d ~.d 
752 " 



7l2 



= 0, 



Vp2 



= 721 = 

Api = —741, 
18.0 GHz, n s0 = v / 274iA s 
and p = 45.9943 Hz. 



B. Incoherent pumping 

In this section, we discuss the possibility to control 
the optically written waveguide structure by means of an 
incoherent pump field. As shown in Fig. EJa), the ab- 
sorption becomes weaker in the two wings than in the 
central area. In practical implementations, this may lead 
to stray fields outside the waveguide region caused by the 
part of the probe field leaking out of the waveguide-like 
structure in the central part into the weakly absorbing 
wings. In order to eliminate this leaking light, we pro- 
pose a spatially dependent incoherent pump field given 
as follows 



p(Ct a 



if (fe) 2 <e- 4 



(12) 
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IV. RESULTS 




FIG. 3. (Color online) Normalized intensity of the probe field 
as a function of propagation distance z. (a) shows the free 
space case, (b) the propagation in a medium prepared with 
a control field and evaluated beyond the paraxial approxima- 
tion, and (c) the corresponding results evaluated in parax- 
ial approximation. The maximum propagation distance is 
zo = 20/^m. The initial width of the probe is w p = 1.0/^m at 
z = 0. Other parameters are as in Fig. [21 



The resulting linear susceptibility is shown in Fig. [2Jb), 
for otherwise same parameters as in (a). It can be seen 
that the central region featuring the waveguide- like struc- 
ture remains the same as in Fig. [2ja) as expected. In 
contrast, in both side wings, where (S7 s /Q s o) 2 < e -4 , the 
incoherent pump vanishes, such that the atoms remain 
in the ground state |1) according to Eq. ©(a). Thus, 
the medium acts as an absorber, and stray-fields outside 
the waveguide region are reduced. In the following, the 
spatially-dependent incoherent pump is used in the nu- 
merical simulations. 



A. Validity of the paraxial approximation 

As a first step, we study two-dimensional (2D) light 
propagation of a Gaussian probe field along z. The probe 
field width is w p = 1.0/mi at z = 0. The medium is pre- 
pared by a Gaussian control laser field f2 s with width 
w s = I/im at z = 0. The respective beam profiles are 
as in Eq. (fTU)) . Note that we do not take into account 
the spreading of the field f2 s . Fig. [3] shows results for the 
normalized intensity of the probe at different positions 
for a propagation distance of zq = 20/im (about 2.59 
Rayleigh lengths Zr ■ = 2itWp/\ p ~ 7.73/im). The three 
panels compare different cases. Panel (a) shows the free- 
space case without control field. Clearly, the diffraction 
spreading of the probe beam can be seen. In (b), results 
are shown for the FDTD simulation beyond the parax- 
ial approximation. It can be seen that the diffraction 
is suppressed by the waveguide structure introduced by 
the control beam. Finally, in (c), results in paraxial ap- 
proximation from the Strang propagation technique are 
shown. While again a waveguide-like propagation with- 
out diffraction, surprisingly, the results differ consider- 
ably from those obtained in (b) beyond paraxial approx- 
imation. In particular, the paraxial approximation leads 
to higher maximum intensity, and the "breathing" -like 
width modulation throughout the propagation is slower 
than in (b). 

Fig. H] shows corresponding sections through Fig. [31 
depicting the intensity profile for the different cases at 
propagation distance zq. In (a), it can be seen that in 
free space, the maximum intensity is considerably atten- 
uated together with a substantial increase in the width to 
w p (zq) ~ 2.85/im. With control field beyond the paraxial 
approximation, the width is slightly narrowed compared 
to the input width to w p (zq) ~ 0.75/im, together with a 
moderate reduction of the maximum intensity. In parax- 
ial approximation, the width is under- and the maximum 
intensity is over-estimated. 

An intuitive way to understand the observed diffrac- 
tionless light propagation is to model the probe as an en- 
semble of rays propagating in a medium rendered by the 
applied control field into a waveguide structure. As the 
probe is propagating in the medium, total reflection at 
the medium boundaries occurs, as the refractive index is 
larger in the central area than in the two wings as shown 
by the solid blue line in Fig. [^a). Thus, the main part 
of the probe intensity is confined in the central area. To 
substantiate this interpretation, we calculated the width 
of the probe as a function of propagation distance z, and 
show the result in Fig.[U[b). It can be seen that the width 
of the probe oscillates periodically throughout the prop- 
agation. These oscillations in the width arise from the 
total refraction. Those rays in the probe which do not 
meet the condition for total internal reflection leave the 
central area and are absorbed in the wings. This leads to 
the reduction in both the energy density and the width 
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FIG. 4. (Color online) (a) Intensity profile of the propagating 
probe field at the medium entrance (dashed black line) and 
output after a propagating distance zo = 20 /im. The dot- 
ted green line shows propagation in free space, the solid red 
line propagation in the medium with control field evaluated in 
paraxial approximation, and the solid blue line corresponding 
results beyond paraxial approximation, (b) Transversal width 
of the probe field as a function of propagation distance z eval- 
uated within (dashed red) and beyond (solid blue) paraxial 
approximation. Parameters are chosen as in Fig. [3] 



of the probe as the propagation distance increases. Note 
that Fig. EJb) again shows the difference in the spatial 
structure of the probe beam with and without the parax- 
ial approximation. 

We thus conclude, that even when confined to a 
waveguide like structure and when propagating only few 
Rayleigh lengths, the paraxial approximation does not 
provide an accurate description of the probe field propa- 
gation. 



B. All-optical switching of light in a branched 
waveguide structure 

We now turn to our main results, on applying the 
waveguide-like structures obtained in the previous sec- 
tion for all-optical switching. For this, we consider light 
propagation in a more complicated branched waveguide 
structure formed by a field fl s consisting of three Gaus- 
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FIG. 5. (Color online) Sketch of the branched waveguide 
structure, (a) shows the spatial profile of Q. s for k = 1. Panels 
(b) and (c) show the corresponding position-dependent real 
and imaginary part of the linear susceptibility in the medium. 
The labels i £ {1, 2, 3} in (a) indicate the three output ports 
corresponding to the i-th control beam Q 3 i. The size of the 
medium is chosen as 60 /im x 200 ftm. The other parameters 
are as in Fig. [2] 



sian laser beams with different geometries defined as 



&s3 = ^soe 



[(g-iEt)< 



(13a) 
(13b) 
(13c) 
(13d) 



The spatial profile of the field ft s is plotted in Fig.[5{a). 
It consists of two parallel beam structures, intersected 
by a tilted one. In this figure, the width w s , the dis- 
placement xq of the parallel beams, as well as the dis- 
placement xt of the tilted beam are chosen as 1.0 /im, 
3.0 /im, and 18.0 /im, respectively, k is a factor by which 



7 



the width of the tilted control beam differs from that 
of the parallel beams. The maximum propagation dis- 
tance is zq = 200 /im, and the tilting angle for f2 S 3 is 
tan# = 2xt/z . Using Eqs. ©, © and (flUj) . we can ob- 
tain the linear susceptibility for the probe, which has real 
and imaginary parts as shown in Fig. [5jb) and [He), re- 
spectively. It can be seen that branched waveguide-like 
structures are generated in the optical response corre- 
sponding to the three beams of fl s . Note that in regions 
in which two of the fields £! s2 and Q s3 overlap, the 
total magnitude of fl s is considerably larger than f2 s0 of 
a single field, which means that the conditions for the de- 
sired maximum reflective index modulation are not sat- 
isfied. This deviation leads to lower spatial dispersion 
with weak gain as shown in the corresponding regions in 
Figs. EJb) and [He). As a consequence, refraction will 
take place when the probe field enters these regions. 

We found that the branched waveguide structure can 
be used to guide the probe light selectively into either of 
the three output ports on the right hand side. For this, 
we consider a Gaussian probe field which is launched into 
the waveguide formed by fi s i left of the intersection point 
with the tilted control beam. As a first step, we have cal- 
culated the light propagation without the tilted control 
field (fi S 3 = 0), such that the probe light remains in the 
upper of the two parallel waveguide structures. We de- 
note this transmitted power after propagation of a single 
optically written waveguide as Pq. We will normalize 
part of the results of the following analysis to this refer- 
ence value Pq, and it is important to remember that Pq 
is smaller than the incident power due to attenuation in 
the waveguide, ions. 

We then switch on the tilted control beam, and cal- 
culate the field power at the three output ports after 
the propagation, which we denote as Pi. From these 
quantities, we calculate the relative transmitted powers 
Ti = Pi/Pa- Here, i = 1 corresponds to the upper paral- 
lel waveguide, i — 2 to the lower parallel waveguide, and 
i = 3 to the tilted waveguide. 

Results are shown in Fig. [6ja-c) for two different tilt- 
ing angles t&n6 £ {0.18,0.26}, as a function of the 
tilted beam width k. Additionally, in (d), the visibil- 
ity V = (Ti - T 2 )/(Ti + T 2 ) is shown, which can be seen 
as a figure of merit for the switching between the two 
parallel output ports. While for k — > 0, all light is emit- 
ted at output port 1 as expected, with increasing k, the 
light is re-routed towards port 2, with a visibility below 
—0.6. This indicates that switching of the output port 
by means of the tilted control field is possible. The corre- 
sponding output probe field structure after propagating 
through the branched waveguide structure without tilted 
field (k = 0) and with k = 4 is shown in Fig. [7] It clearly 
shows that the probe pulse in output 2 exceeds that in 
port 1. 

Fig. [7] also shows that the output power can be larger 
than Pq. This is also shown in Figs.[5ta-c). The origin of 
this increase is the slight gain in the field overlap regions. 
It should be noted, however, that overall the sum of the 





FIG. 6. (Color online) (a,b,c) Transmitted powers Ti, T2, and 
T3 at the three output ports indicated in Fig. [5ja) after the 
propagation through the branched waveguide medium, (d) 
shows the visibility V = (Ti -T 2 )/(Ti +T 2 ) for the switching 
of light between the two output ports 1 and 2. All results are 
plotted against the relative width of the tilted control beam 
k. The blue solid (red dashed) line shows results for tilting 
angle tan 9 = 0.18 (tan (9 = 0.26). The probe beam initially 
has a Gaussian shape with w p = 1.0 /J.m. Other parameters 
are as in Fig. [2] 



total transmitted probe intensity in all three output ports 
is lower than the input power. Values larger than Po 
only indicate that the transmitted power is larger than 
the power transmitted through a single optically written 
waveguide, as part of the absorption is compensated by 
the gain in the overlap region. 

In the following, we provide a simple explanation for 
the observed light switching to different output ports. If 
the tilted beam is switched off, the probe field propagates 
without diffraction to output port 1 as discussed in the 
previous section. When the tilted beam S 3 is turned 
on, and if suitable parameters are chosen for the tilting 
angle, the tilted beam width, and the atomic density N, 
the probe field can be refracted into the tilted beam f2 S 3. 
To understand this, on should note that the region in 
which two of the three control fields overlap has a real 
part of the susceptibility varying with the propagation 
direction z. Due to this variation, as well as due to the 
missing guiding in perpendicular direction in the over- 
lap region, parts of the probe field can be redirected into 
the tilted beam. Starting from a tilted beam of vanish- 
ing thickness, initially the proportion of redirected light 
increases with increasing thickness. After the first inter- 
section point, at the second intersection point, in turn 
light can be redirected out of the tilted beam into the 
lower parallel control beam. 

Note that the transmitted power T\ at output port 1 
oscillates with the thickness of the tilted control beam 
determined by k. This is likely an interference effect. 
As discussed above, the refractive index varies along the 
propagation direction z in the waveguide formed by Q s \ 
in the intersection point. When the probe field enters 
this area, some part of it is refracted back and forth 
repeatedly along the z direction. Depending on k, the 
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FIG. 7. (Color online) Spatial field configuration in the plane 
transverse to the propagation direction at the output of the 
medium at zo = 200 fim. The red solid line shows results 
without tilted control beam, whereas the blue dashed line 
shows results including a tilted control with relative width k = 
4. The vertical green dashed lines indicate the positions of the 
three output ports. The tilting angle is chosen as tan# = 0.18. 
Other parameters are as in Fig. 
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FIG. 8. (Color online) Visibility V = (Ti - T 2 )/(T 1 + T 2 ) 
as a function of the tilting angle of the control field for three 
different relative widths of the tilted beam. The upper red 
line shows k = 2, the middle blue line k = 3, and the lower 
green line k = 4. Other parameters are as in Fig. [6] 

interference of the different channels possible for the 
light changes between (partly) constructive or destruc- 
tive, thus resulting in the oscillations in T\. 

Qualitatively, one would expect that the probe field 
can more easily be redirected from H a i into the tilted 
beam fi s3 and subsequently into the lower parallel guide 
51 S 2 if the tilting angle 8 is smaller, and if the atomic 
density N is larger. We have verified this by plotting the 
transmitted power Tj and visibility V versus the width of 
fi S 3 also for a larger tilting angle tan 9 = 0.26 in Fig. [5] 
It can be seen that the probe transferred into the two 
output ports 2 and 3 is generally smaller at tan0 = 0.26 
as compared to that at tan0 = 0.18. We then further 
calculated the visibility V as a function of the tilting an- 
gle tan (9 for several different widths of the tilted control, 
results are shown in Fig. [5] It can be seen that V in- 



FIG. 9. (Color online) The effect of the incoherent pump on 
the output of the probe. The angle and width of the probe 
are chosen as tan# = 0.18 and k = 4 respectively. Other 
parameters are as in Fig. [6] 

creases as the tilting angle 9 increases for all widths of 
the tilted beam. 

It should be noted that since the tilting angle in our 
case is relatively large (tan 9 £ {0.18, 0.26}), the paraxial 
approximation cannot be applied. This is clearly indi- 
cated by the non-vanishing transmitted power T3 which 
arises due to effects beyond the paraxial approximation 
described by the term (ic/2u} p )d 2 E p /dz 2 in Eq. (jHJ). 

As discussed in Sec. II V A} since we have chosen a rela- 
tively weak incoherent pump, the probe will be absorbed 
as it propagates in the medium. This is an obstacle in 
experimentally realizing the proposed scheme. However, 
the output of the probe can be improved by applying 
stronger incoherent pumping. Fig. |9] shows the transmit- 
ted power Ti(i = 1,2,3) and visibility V as a function 
of the incoherent pump rate p. It can be seen that the 
transmitted powers Tj can be sensitively controlled via 
the intermittent gain introduced via p over a wide range 
of output powers. Note that due to the slightly longer 
propagation distance to output port 2 compared to the 
path to 1, the transmitted power T2 grows faster than T\ 
as the incoherent pump increases. For the same reason, 
T3 grows even more rapidly. It is important to note that 
while the transmitted powers can be controlled, the visi- 
bility for the switching remains approximately the same 
over the whole range of p. This suggests that the visi- 
bility is mainly determined by the geometry of the op- 
tically written structure, consistent with the waveguide 
interpretation. In contrast to the absorption and gain 
properties, this geometry is largely independent of the 
pumping p, with the exception of slight changes in the 
beam profiles with p. 



V. CONCLUSION 

We have investigated light propagation in an opti- 
cally written waveguide structure. The waveguide struc- 
ture is prepared in a medium tailored with coherent 
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control fields and an incoherent pump field. A spa- 
tially dependent Gaussian control field is applied to cre- 
ate a waveguide-like spatial refractive index variation 
within the medium, such that a probe beam can prop- 
agate within the waveguide structure essentially with- 
out diffraction. Our initial calculations show that al- 
ready in a single optically written waveguide, accurate re- 
sults for the beam propagation cannot be obtained within 
the paraxial approximation. To analyze the propagation 
dynamics, we therefore numerically integrate the corre- 
sponding equations without paraxial approximation for 
the probe fields. 

As our main result, we have demonstrated that a con- 
trollable branched-waveguide structure can be optically 
formed inside the medium by applying a spatially depen- 
dent control field consisting of two parallel and one tilted 
Gaussian beams. The tilted beam can be used to selec- 
tively steer a probe beam between two different output 
ports. In the overlapping regions of the Gaussian control 
beams a rapidly varying refractive index is created that 



guides the probe beam to propagate in a particular out- 
put branch. Increasing the amplitude of the incoherent 
pump beam can be used to reduce the absorption of the 
medium. 
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Appendix A: Appendix 

The explicit expressions for A,B in Eqs. © are as 
follows: 



A = -i (-732 - i (A c - A p2 - A s )) ((-*7sa + A c - A p2 - A s ) (i l52 + A p2 + A s ) + il 2 c - flf) (743 + i (A c - A p2 ) + 

0| &s ((*74 3 - A c + A p2 ) (-z 732 + A e - A p2 - A s ) - + Q 2 S ) \ 

732 + i (Ac - A p2 - A s ) (732 + i (A c — A p2 - A.)) ((-vyaa + K - A p2 - A.) (£752 + A p2 + A.) + ftjj _ Q2) j 

(Ala) 

B = ( 732 + 7 54 + i (A c - 2A p2 - A,)) Q 2 S ((i 743 - A c + A p2 ) (-i 732 + A c - A p2 - A s ) - fl 2 c + Q 2 S ) + 
((-2732 + A c - A p2 - A s ) ( 2752 + A p2 + A.) + Q 2 C - Q 2 S ) ((732 + * (A c - A p2 - A,)) x 

((-*743 + A c - A p2 ) (i 754 + A p2 ) + n*) + (754 - iA p2 ) n 2 s ) (Alb) 
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